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The properties of narrow- line Seyfert 1 (NLS1) galaxies, the links and correla- 
tions between them, and the physics behind them, are still not well understood. 
Apart from accretion rates and black hole masses, density and outflows were 
speculated to be among the main drivers of the NLS1 phenomenon. Here, we 
utilize the diagnostic power of the [S II] AA6716, 6731 intensity ratio to measure 
the density of the narrow-line region (NLR) systematically and homogeneously 
for a large sample of NLS1 galaxies, and we perform a comparison with a sam- 
ple of broad-line typel Active Galactic Nuclei (AGN). We report the discov- 
ery of a ' 'zone of avoidance 1 in density in the sense that AGN with broad lines 
(FWHM(H / 9) > 2000 km s" 1 ) avoid low densities, while NLS1 galaxies show a 
wider distribution in the NLR density, including a significant number of objects 
with low densities. A correlation analysis further shows that the Eddington ratio 
L/LEdd anti-correlates with density. We investigate a number of different models 
for the 'zone of avoidance' in density. Supersolar metallicities and temperature 
effects, a strong starburst contribution in NLS1 galaxies, different NLR extents 
and selective obscuration are considered unlikely. Possible differences in the frac- 
tion of matter-bounded clouds and differences in the interstellar media of the 
host galaxies of NLS1 galaxies and broad-line Seyfert 1 (BLS1) galaxies can only 
be tested further with future observations. We tentatively favor the effects of 
winds/outflows, stronger in NLS1 galaxies than in BLS1 galaxies, to explain the 
observations. 

Subject headings: galaxies: density - galaxies: ISM - galaxies: emission lines - 
galaxies: active - galaxies: Seyfert 

1. Introduction 

Optical and X-ray observations over the past few decades revealed a new sub-class 
of AGN, termed Narrow-line Seyfert 1 (NLS1) galaxies (e.g., Gaskell 1984; Osterbrock & 
Pogge 1985). NLS1 galaxies are intriguing due to their extreme emission line and continuum 
properties. Their optical broad lines are narrower (FWHM FL3 < 2000 km s -1 ) than in 
'normal' broad-line Seyfert 1 (BLS1) galaxies and they show strong Fe II emission. Their X- 
ray spectra are sometimes, but not always, very soft (e.g., Zhou et al. 2006, and references 
therein). Many emission-line and continuum properties of AGN were found to correlate 
strongly with each other (e.g., Boroson & Green 1992, BG92 hereafter; Wang et al. 1996; 
Lawrence et al. 1997; Grupe et al. 1999; Vaughan et al. 2001; Xu et al. 2003; Sulentic et al. 
2000, 2003; Grupe 2004). The strongest variance, often referred to as 'Eigenvector!' (EV1), 
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is defined by the correlation between the width of the H/3 emission line and the strength 
of the [O IIl]/H/3 emission line ratio, and the ant i- correlation with the Fe II/H/5 ratio (e.g. 
BG92). NLS1 galaxies are placed at one extreme end of EV1 parameter space. The most 
common interpretation is that this regime is goverened by the highest Eddington accretion 
rates and/or lowest black hole masses (e.g., BG92; Sulentic et al. 2000). 

Among other parameters^, the density of an outflowing wind was firstly speculated to 
be a prominent driver of EV1 by Lawrence et al. (1997), given the connection of Fe II 
strength with the presence of low-ionization, blueshifted broad absorption lines, and with 
blue- asymmetric emission lines. Winds and outflows play a crucial role in understanding the 
physics and evolution of AGN (e.g., Elvis 2000, Hopkins et al. 2005), and there is ample 
observational evidence for winds and outflows in AGN from sub-kpc to galactic scale (see, 
e.g., Sulentic, Marziani & Dultzin-Hacyan 2000; Veilleux, Cecil, & Bland-Hawthorn 2005, for 
recent reviews). In enriching the nuclear environment with matter from the central region, 
winds may have an important impact on the gas densities in the emission-line regions. 

Regarding NLS1 galaxies, their high ratios of L/Lgdd are likely particularly efficient 
in driving outflows. The large [O III] A5007 blueshifts observed in some NLS1 galaxies 
are interpreted straightforwardly as the result of an out flow (e.g., Zamanov et al. 2002; 
Aoki et al. 2005, Boroson 2005). iDewangan et al.l (120011 ) suggested that stronger outflows 
would likely push the broad-line region (BLR) further radially outward thereby resulting in 
narrower H/3 lines in NLS1 galaxies. They further speculated that the observed flux ratios 
from the NLR and BLR can be explained in terms of density enhancements (see also Wills 
et al. 2000). Outflows in NLS1 galaxies have been observed in terms of both blueshifted 
UV absorption lines (e.g., Laor et al. 1997a,b; Goodrich 2000) and UV emission lines (e.g., 
Leighly & Moore 2004). 

Regarding the density of the emission-line regions of NLS1 galaxies, only few previous es- 
timates exist (e.g., Wills et al. 2000; Kuraszkiewicz et al. 2000; Sulentic, Marziani & Dultzin- 
Hacyan 2000; Marziani et al. 2001; Bachev et al. 2004). Thos e which do exist, actually 



lead to partially conflicting results. While some sample studies (jKuraszkiewicz et all 12000 



1 See lKomossa et al] (|2006b| ) for a recent discussion. Briefly speaking, other parameters considered to be, 
perhaps, relevant in explaining NLSf properties are orientation (e.g., Osterbrock & Pogge 1985; Puchnarewicz 
ct al. 1992; Collin et al. 2006; Zhang & Wang 2006), the physics behind NLSl's radio properties (e.g. 
Komossa et al. 2006b), metallicity (e.g., Mathur 2000; Komossa & Mathur 2001; Nagao et al. 2002; 
Shemmer & Netzer 2002; Romano et al. 2004; Fields et al. 2005), (ionized) absorption (e.g., Komossa & 
Meerschweinchen 2000; Gierlinski & Done 2004), and the galaxies' location on the Mbh — cr plane (e.g. 
Mathur et al. 2001; Wang & Lu 2001; Grupe & Mathur 2004; Bottc ct al. 2005; Mathur & Grupe 2005; 
Komossa & Xu 2007) 
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Wills et alJbOQOl ; iMarziani et al.ll200ll ; iBachev et al.ll2004j ) suggested a high-density BLR in 
NLS1 galaxies, based on the large Si III] A1892/C III] A1909 ratios measured from UV spectra 
@, other inquiries favored low density emission-line clouds. Rodriguez- Ardila et al. ( 2000a ) 
studied the emission-line properties of a sample of Seyfert 1 galaxies, including 7 NLS1 galax- 
ies using optical and near-IR spectroscopy. They found that the ty pical density of the [S II 



emitt ing zone is lower in NLS1 galaxies than in BLS1 galaxies. iRodriguez-Pascual et al. 



(119971 ) tentatively favored a low-density BLR i n NLS1 galaxies, inferred from their pho- 
toionization modeling of the UV emission lines. iFerland fe Perssonl (119891 ) derived densities 
similar to the canonical BLR value in st rong Fe II emit t ers, b ased on the ratio of the for- 
bidden lines relative to the Ca II triplet. iBaskin fe Laorl (120051 ) reported that the density of 
the [O III] A5007 emitting gas in the NLR decreases with steeper observed soft X-ray slope. 



In addition, a few individual objects were inspected more closely (ILaor et al.l 11997b 
Leighlv fe Moore! 120041 ; IVeron-Cettv. Jolv fe VerorJ 12004! ; fveron-Cetty et alJbood ). In par- 



ticular, analyses of the high S/N UV spectrum of the prototype NLS1 galaxy, IZwl, ind i- 
cated a high BLR density o f 10 11 c m" 3 and NLR density of 5 x 10 5 cm -3 (ILaor et al.lll997bl ). 



Veron-Cetty. Joly fe Veronl (120041 ) found on the same object from optical spectroscopy that 



the density of the low- ionization part of the NLR is of the order of 10 6 7 cm 3 . 

Given that few systematic measurements exist at all, and that those which do produce 
partially conflicting results, it is important to explore this topic further. We present for the 
first time a study of the NLR density for one of the largest homogeneously analyzed NLS1 
samples to date and compare it with that of BLS1 galaxies. The electron density of the 
NLR can be measured by making use of the density-sensitive line ratios [O II] A3729/A3726 
and [S II] A6716/A6731 (e.g., Osterbrock 1989). In practice, the most important density 
diagnostic is the [S II] line ratio. Usually the density can not be directly inferred from the 
[O II] line ratio because the [O II] AA3726, 3729 doublet is often unresolved. More indirectly, 
other line ratios will also change with density. For instance, a higher-density NLR would 
strongly boost the [O I] A6300 line, thus the intensity of [O I] can be used to probe higher 
density regions (e.g., Komossa & Schulz 1997; Barth et al. 2001). 

This work is part of a series of papers investigating the properties of the emission-line 
regions of type 1 AGN, including NLS1 galaxies and BLS1 galaxies. In this first paper of the 
series, we focus on topics related to density. Specifically, we attempt to answer the following 
key questions: (1) is there any difference in the NLR density between NLS1 galaxies and 
BLS1 galaxies? (2) If so, do trends in density correlate with other parameters? and (3) what 



2 One of the a lternative explanation s is that the carbon has been removed from gas phase by depletion 
onto dust grains |Crenshaw et al.ll2002h . 
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are the key physical drivers to explain differences in the NLR of NLS1 galaxies and BLS1 
galaxies? 

This paper is organized as follows. We present the data base and the sample selection 
in Sect. 2. In Sect. 3, we describe our method of the optical spectral analysis. The sample 
classification and an investigation of selection effects is provided in Sect. 4. The key result, 
the detection of a 'zone of avoidance 1 in the NLR density, is reported in Sect. 5. In Sect. 6 we 
discuss the reality of the zone of avoidance, followed by a discussion on its origin in Sect. 7. 
We summarize our conclusions in Sect. 8. 

We use the terms NLS1 galaxies and BLS1 galaxies collectively for high-luminosity and 
low-luminosity objects, i.e., NLS1 galaxies for narrow-line type 1 quasars and for narrow-line 
type 1 Seyfert galaxies, and BLS1 galaxies for broad-line type 1 quasars and for broad-line 
type 1 Seyfert galaxies, respectively. Throughout this paper, a cosmology with H = 70 km 
s -1 Mpc" 1 , f^M = 0.3 and Qa = 0.7 is adopted. 



2. Sample selection 

The combination of X-ray and optical observations has proven to be an efficient way in 
detecting NLS1 galaxies over the last decades (e.g., see Pogge 2000 and Veron-Cetty, Veron 
& Gongalves 2001, hereafter VVG01, for reviews). The uniform optical spectroscopic galaxy 
survey known as the Sloan Digital Sky Survey (SDSS, York et al. 2000) is an excellent 
data base which enables us to study AGN properties and, in particular, optical emission line 
properties in a homogeneous way. 

In this paper, in order to measure NLR densities, we homogeneously analyze and com- 
pare optical emission line properties of a large number of NLS1 galaxies with BLS1 galaxies 
based on spectra obta ined in the course of the SDSS. We use the 3rd data release, DR3 



(lAbazajian et al.ll2005l ). The DR3 spectroscopic program covers an area of about 4188 deg 2 . 
The two double spectrographs produce data covering a wavelength range 3800-9200A at 
a spectral resolution rj 2000. Exposure times for spectroscopy are determined such that a 
signal-to-noise (S/N) of at least 4 pixel -1 at g = 20.2 is reached. The SDSS spectroscopic 
pipeline (Stoughton et al. 2002) performs spectral extraction, sky subtraction, removal of 
the atmospheric absorption bands, wavelength and flux calibration, and estimates the error 
spectrum. The processed DR3 spectra have not been corrected for Galactic extinction, but 
spectrophotometric calibration has been considerably improved since the First Data Release 



(DR1). We refer the reader to lAbazajian et al.l (120051 ) for details of the changes 
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2.1. The NLS1 sample 

First, we extracted all NLS1 galaxies included in the the 11th edition of the "Catalogue 
of Quasars and AGN" (Veron-Cetty & Veron 2003; VV03 hereafter.). The defining criterion 
of NLS1 galaxies in VV03 is FWHM(H/3) < 2000 km s~ x . We further take into account only 
NLS1 galaxies with redshift z less than 0.3. The redshift cut is imposed to ensure that [S II] 
AA6716,6731 is observable, and furthermore, to ensure that the doublet is in a region free of 
noise from strong night-sky emission lines. The list of 309 selected NLS1 galaxies was then 
cross-correlated with the SDSS DR3, in order to obtain a homogeneous set of spectra of a 
large sample of NLS1 galaxies for uniform spectral analyses. This procedure resulted in the 
selection of 119 sources. In order to get an accurate density measurement, we require [S II] 
AA6716,6731 to have S/N greater than 5. We also removed a handful of objects with spectral 
defects (e.g. problems with sky subtraction, or missing signal over a range of wavelengths; 
see Strateva et al. 2003) by visual inspection. This leaves a total of 58 objects. 

2.2. The BLS1 sample 

Recently, Boroson (2003) presented a sample of 107 low-redshift (redshift z < 0.5) 
type 1 AGN from the SDSS Early Data Release (EDR, Stoughton et al. 2002)|. 82 out of 
these objects match our redshift constraint, i.e., z < 0.3. We use this sample to build the 
control BLS1 sample imposing the same redshift cut and S/N limit of [S II] AA6716,6731 (i.e., 
z < 0.3 and S/N > 5) as to the NLS1 sample. The processed DR3 spectra were downloaded 
for spectral inspection in accordance with the data set of the NLS1 sample. With these 
restrictions, 48 objects survive, 13 of which overlap with the NLS1 sample. 

3. Spectral Analysis 

The optical emission-line properties (line widths and line ratios) are among the defin- 
ing criteria of the NLS1 phenomenon. The exact and homogeneous measurements of the 
emission-line parameters of AGN provide us with basic knowledge about emission-line re- 
gions, and enable us to investigate the relationship between NLS1 galaxies and BLS1 galaxies. 
However, the published data were compiled in a heterogeneous way. The use of different line 



3 We do not div i de th e Seyfert 1 class into subclasses such as Seyfert 1.5-1.9 (as defined by 
Osterbrock fc Poggej (|1985h ). since the de finitions of Seyfert classifications d epend on the resolution of the 
spectra used and the noise in the spectra ( Grupe et al. 19991 : Goodrich 1989f) . 
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profiles (e.g., Lorentzian vs Gaussian profiles) leads to strong differences in the H/3 widths 
(e.g., VVG01). In addition to the Lorentzian and Gaussian representations of line profiles, 
a direct measurement of the width at half of their maximum intensity is sometimes adopted 
(e.g., Williams, Pogge & Mathur 2002; Boroson et al. 2003). Also, the strong Fe II contam- 
ination makes it difficult to measure H/3 and brings large uncertainties in determining the 
H/9 line width. 

Here we perform a homogeneous spectral analysis for the NLS1 and BLS1 samples. 
The objects are then re-grouped into the NLS1 or BLS1 sample for further investigation, 
according to the widths of the broad H(3 line (Sect. 3.2) obtained in our spectral analysis. 

In a first step, the SDSS spectra were corrected for the Galactic extinction using the 
reddening map of Schlegel et al. (1998), and then shifted to their rest wavelength, adopting 
the SDSS value of the redshift from the header of each spectrum. Once these steps were 
completed, we removed the stellar continuum, subtracted the Fe II complexes and performed 
the spectral analysis for the emission lines following the procedures described below. 

3.1. Subtraction of starlight and nuclear continuum 

SDSS spectra are acquired with a pair of fiber-fed spectrographs. Each fiber subtends a 
diameter of 3", corresponding to ~ 6.5 kpc at z — 0.1. This aperture is large enough to let 
through not only the emission from the nucleus, but also a substantial amount of starlight 
from the host galaxy (e.g., Vanden Berk et al. 2001, Hao et al. 2005). The accurate removal 
of the stellar contribution is essential to reliably measure the emission-line spectrum, such as 
the line widths and line strengths. Particularly, the reliable classification of NLS1 galaxies 
is strongly dependent on the measurement accuracy of the width of the broad H/5 line. 
Furthermore, in many of the spectra there is a clear contribution from blends of Fe II line 
emission. Well-studied NLS1 galaxies usually show strong optical Fe II emission features on 
both the blue and red sides of the H/3-[0 III] complex (e.g., Veron-Cetty & Veron 2003). 
In order to reliably measure line parameters, we choose those wavelength ranges as pseudo- 
continuum, which are not affected by prominent emission lines, and then decompose the 
spectra into the following 4 components (see Zhou et al. 2006 for details): 

• A starlight component modeled by 6 synthesized galaxy templates, which were built 
from the synthetic spectral library of Bruzual & Chariot (2003) using the algorithm 
of Ensemble Learning for Independent Component Analysis (EL-ICA, Lu et al. 2006). 
These templates were broadened by convolution with a Gaussian to match the stellar 
velocity dispersion of the host galaxy. 
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• A power-law continuum to describe the emission from the active nucleus. 

• An Fe II template obtained by Veron-Cetty & Veron (2003). This template covers the 
wavelengths between 3535—7534 A, extending further to both the blue and red wave- 
length ranges than the Fe II template used in BG92. This makes it more advantageous 
in modeling the Fe II emission in the SDSS spectra. We assume that Fe II has the 
same profile as the broad component of H/3 (see the next subsection). 

• A Balmer continuum generated in the same way as Dietrich et al. (2003). 

The modeling is performed by minimizing the reduced y 2 in the fitting process. The 
final multi-component fit is then subtracted from the observed spectrum. An example of the 
residual spectrum is plotted in Fig. 1. 



3.2. Decomposition techniques 



The multicomponent-subtracted spectra are used to measure the non-Fe II line proper- 



ties. The broad Balmer lines in AGPs 



Miller et al 



metric ( 



range in width (IQsterbrock fe Shuder 



1992; 



exhibit a wide variety of profile shape s and a large 



1982; 



de Robertislll985l : ICrenshawl[l9~86l : Istirpei Il991 



Veron-Cetty. Veron fe Goncalvesl 2001 ). and they are often strongly asym 



Corbinl Il995l ) . In many cases the Balmer lines are mixtures of broad and narrow 



components. Differences in the rela tive strengths of th e se components account for much of 



the diversity of br o ad line profiles (jFrancis et~aL 



1992 



Wi 



1994 ; ICorbinlll995l . Il997l : IVeron-Cettv. Veron fe Goncalves 



Is et al. 1993; Brotherton et al. 



200 ll ). A proper decomposition 



of the NLR and BLR line emission contribution of the Balmer lines is of great importance 
to address the physical properties of the emission-line regions. Particularly, the width of the 
broad Hf3 profile has significant impact on the reliable classification of a galaxy as a NLS1 
galaxy. 

In order to measure the parameters of the BLR emission lines, the NLR line emission 
contribution has to be removed first. Using homogeneous sets of spectra, previous studies 
(e.g., Filippenko & Sargent 1988; Ho et al. 1997; Greene & Ho 2004; Zhou et al. 2006) have 
shown that the narrow Balmer emission profiles are well matched to those of [S II] AA6716, 
6731. Moreover, the widths of [S II] and [N II] doublets trace the stellar velocity dispersion 
of galaxy bulges better than that of [O III] within the uncertainties (Greene & Ho 2004). 
Therefore, we use the strong profiles of [S II] AA6716, 6731 (those with S/N >5) as a NLR 
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template profile for the narrow component of the Balmer linefl 



To characterize the NLR emission-line profiles, we fit these lines using Gaussian profiles. 
Most [S II] AA6716,6731 and [N II] AA6548,6583 lines can be well fit employing one single 
Gaussian profild^]. A large fraction of the [O III] AA4959,5007 lines show strongly asymmetric 
profiles. In those second component is then added to represent the line wings. When 

we fit [S II] AA6716,6731, [N II] AA6548,6583 and [O III] AA4959,5007 lines, the separation 
of the lines of each doublet is fixed to the laboratory value. Each pair is assumed to have 
the same profile^]. While the flux ratios of the [N II] and [O III] doublets are fixed to the 
theoretical ratio of 3 : 1, the intensity ratio of the [S II] doublet is measured, and then used 
to derive the density. 



We used the IRAF package SPECFIT (jKrissJ |l994j) to measure blended lines and sep- 
arate the NLR from the BLR emission. The actual profile shape of the broad Balmer lines 
of NLS1 galaxies is still an issue debated. While the broad Balmer line component can be 
well represen ted by Gaussian profiles, particularly by a combination of multiple Gaussian 
compo ne nts ( Rodriguez- Ardila et al. 2000bl; Nagao et al. 2002 ; Xu et al. 2003 ; Greene fc Ho 
2005al lbl; iDietrich. Grenshaw fc Kraemerl 2005^ , successful fitting can also be accompli shed 



with a single Lorentzi an profile (ILeighlv 



1999 



Veron-Cettv. Veron fc Goncalved 1200 ll ) for 



some NLS1 galaxies. IVeron-Cetty. Veron fc Goncalved (120011 ) suggested that the broad 



Balmer lines of NLS1 galaxies were better fitted with a single Lorentzian profile than a 
single Gaussian profile (see also Sulentic et al. 2002), since many NLS1 galaxies were lo- 
cated in the H II regio n in th e ir dia gnostic diagrams if a single Gaussian profile was adopted. 
However, as noted in lEvand (119881 ). the choice of Gaussian or Lorentzian profiles as repre- 
sentatives of the observed emission lines may bear no physical meaning. Furthermore, the 
broad-line profiles show complex shapes and asymmetries that can not be described with a 
single component, indicating t he presence of two or multiple compo nents, independent from 
the type of profile used to fit (IDietrich. Grenshaw fc Kraemerl 120051 ). 



In order to isolate the narrow and broad components of the Balmer lines, we fix the 



4 Occasionally, we also make use of spectra with weak or absent [S II] emission. These are used to measure 
[O I] (see Sect. 5.3). In such a case, the [O III] profile is used as a substitute for [S II]. 

5 Using spatially-resolved HST spectra, Rice et al. (2006) identified both blue and red asymmetries in 
some [S II] line profiles, which are primarily due to nuclear line-emitting gas, rather than more symmetric 
emission from the NLR on larger scales. 

6 The line width of each component of the doublet may be slightly different due to stratification in the 
NLR. However, Rice et al. (2006) found the widths of the two lines differ by less than 3%, which is within 
the error of the Gaussian fit parameters and hard to discern with data of moderate S/N. 
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width of the narrow component (determined from the width of [S II] as described above) 
and only leave its strength as fit parameter. The broad part of the line profiles we fit by 
using a combination of two Gaussian profiles, as well as a single Lorentzian profile. For the 
approach with a multi-Gaussian components fit, we measure the FWHM of the profile and 
its integrated flux from the final combined model for the broad component. The individual 
Gaussian components have no physical significance by themselves but are only used to serve 
as a description of the complex line shapes of the broad components as far as the data quality 
allows. While comparable and equally reasonable results can be achieved with both single 
Lorentzian and two Gaussian profiles for most broad components of the Balmer lines of NLS1 
galaxies, generally no acceptable fit is possible when employing the Lorentzian profiles to fit 
BLS1 galaxies. We will use the results of the multiple Gaussian fit for NLS1 galaxies for 
classification and further correlation analysis, for its simplicity, and furthermore, for a direct 
comparison with the BLS1 control sample and with previous studies (e.g. BG92; Grupe et 
al. 1999; Vaughan et al. 2001). However, for comparison purposes, we still report results 
from Lorentzian fitting in our key figures (Fig. 5, 6). We further checked the location of 
NLS1 galaxies and BLS1 galaxies in the diagnostic diagrams using the standard emission-line 
ratios. In contrast to Fig. 8 of VVG01, most of our NLS1 galaxies locate in the AGN regime 
for both Gaussian and Lorentzian fits. 

The broad component of H/5 is referred as H/3 b , and the narrow component as H/3 n . The 
dominant uncertainties in the line parameter measurements can result from the continuum 
subtraction and component decomposition. The latter depends on the line profile shape. 
Although the line measurements of Gaussian and Lorentzian fits are precise, they are only 
true estimates if the lines can be correctly represented by the profiles. The average error 
in the width of narrow lines is of approximately 5%. Errors in the width of broad lines, 
introduced by the different profile types, e.g., two Gaussians or a single Lorentzian, are 
about 10%. For most objects, the uncertainties of flux measurements of emission lines are 
less than 10%. The typical measurement error of the [S II] doublet is about 5% for our 
sample. 



4. Sample re-classification and selection effects 

4.1. Re-classification of the sample 

In order to follow the standard practice to distinguish between NLS1 galaxies and BLS1 
galaxies by FWHM(H/3 b ), we re-classified the objects of our sample accordingly, after having 
carefully measured FWHM(rL3 b ). A total of 55 objects with FWHM(H/? b ) < 2000 km s _1 are 
included in the final NLS1 sample, while 39 with FWHM(H/3 b ) > 2000 km s" 1 are included 
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in the BLS1 sample. 

4.2. Luminosity and redshift distributions 

As described in Sect. 2, we do not set any selection criteria upon our samples, which 
combines data from various different sources and thus is not statistically complete. It is 
therefore necessary to examine whether or not the NLSI and BLS1 samples are drawn from 
comparative populations such that luminosity effects (in particular, the number of quasars 
vs Seyferts in each sample) and evolutionary effects can be excluded. In order to check the 
possible biases between the two samples, we look into the absolute magnitude distribution 
and the redshift distribution. 

First, we investigate whether or not there is a systematic difference in the luminos- 
ity distribution between the two samples. Here we trace the luminosity with the absolute 
% magnitude, which is calculated from the i_psf magnitude (e.g. Schneider et al. 2005; Van- 
den Berk et al. 2004), by correcting the % measurements for Galactic extinction (Schlegel 
et al. 1998) and assuming a powerlaw spectral energy distribution (SED) (f v oc u a ), where 
a = —0.5 (e.g., Vanden Berk et al. 2001). Following the previous studies of the SDSS AGN 
(Schneider et al. 2003, 2005; Nagao et al. 2005), we define objects that have luminosities 
larger than Mj = —22.0 as quasars. We show the histograms of the absolute i magnitudes 
for the NLSI galaxies and BLS1 galaxies in Fig. 3. We apply the Kolmogorov-Smirnov (K-S) 
statistical test to check how closely the distributions resemble each other. The resultant K-S 
test probability is 0.33, which means the two distributions are statistically indistinguishable. 

Secondly, we explore the redshift distribution. The histograms of the sources' redshifts 
are shown in Fig. 3. The NLSI and BLS1 samples possess similar average redshifts. The 
value and la deviations for the two samples are 0.128 ±0.068 and 0.142 ±0.066, respectively. 
Moreover, since both samples are compiled from low-redshift objects with z < 0.3, no strong 
redshift bias is expected, anyway. We conclude that our two samples do not show significantly 
different luminosity and redshift distributions. 

4.3. Optical Fe II emission versus [O III] emission 

In this section we wish to demonstrate that our (NL)S1 sample follows the same cor- 
relations found previously for (NL)Sls as a class, and therefore, that our NLSI galaxies are 
representative for the NLSI population as a whole, as far as previously known correlations 
are concerned. Among optical emission lines of NLSI galaxies and BLS1 galaxies, the most 
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striking correlation is the so called EV1 (e.g., BG92; Grupe et al. 1999; Sulentic et al. 2002, 
2003; Xu et al. 2003), i.e., that [O III]/H/3 emission is weak in objects with strong Fe II 
emission and vice versa. NLS1 galaxies are at the extreme end of this correlation, often 
showing strong Fe II emission and weak [O IIl]/H/3 tota i emission. 

The optical Fe II strength, R4570, is defined as the ratio of the flux of the Fe II complex 
between the rest wavelength 4434A and 4684A to the total H/3 flux, including the narrow 
component (e.g. BG92, VVG01). The average values of R4570 for our NLS1 and BLS1 
samples are 0.75 and 0.30, respectively. The regime of high values of R4570 is solely occupied 
by NLS1 galaxies. The results are consistent with previous studies (e.g., Rice et al. 2006 and 
references therein). In Fig. 4, we plot R5007, the ratio of [O III]A5007 flux to the total B./3 
flux, as a function of R4570. Both emission lines span almost the whole range from the low to 
high end observed in NLS1 galaxies and BLS1 galaxies (e.g., Fig. 11 in VVG01), suggesting 
there is no preferential selection effect to objects with extreme properties. Moreover, an anti- 
correlation between the two parameters is present, with a correlation coefficient r s = —0.37. 
The probability of the null correlation is P nu n < 10~ 3 . These quantities are calculated from 
the Spearman rank-order correlation analysis. Hence, the NLS1 and BLS1 samples under 
study can be considered as representative samples of these two subclasses, though they are 
not statistically complete. 



5. A zone of avoidance in the NLR density 
5.1. Measurement of density 



We use the density diagnostic [S II] A6716/ A6731 to m e asure the NLR electron den- 
sity. The electron density is ca lculated following lOsterbrockl (119891 ) using the task temden 
of the IRAF package STSDAS (IShaw fe Dufourlll994j ). This approach comes with two lim- 
itations. Firstly, around densities of 10 3 -10 4 cm -3 we approach the critical densities of 
[S II] A6716 (1.5 x 10 3 cm" 3 ) and [S II] A6731 (3.9 x 10 3 cm- 3 ). Secondly, the [S II] line ratio 
also depends on temperature. However, the dependence is only weak in the range of tem- 
peratures of the regions studied. We fix the electron temperature at T c = 10, 000K, typical 
for photoionized gas in the NLR, and comment on possible temperature effects later. The 
[S II]-based density measurement allows us to probe the [S Il]-emitting part of the NLR. In 
order to get a first impression also on the high-density part of the NLR, later we make use 
of the strength of the [O I] line (Sect. 5.3)0. 



7 We note that the NLR is most likely composed of emission-line clouds with a range in densities (e.g., 
Komossa & Schulz 1997; Peterson 1997; Rodriguez- Ardila, Pastoriza & Donzelli 2000a; Brinkman et al. 
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The [SII] ratio is only a good density diagnostic in a certain range of densities because 
depopulation of the upper levels changes from being primarily radiative at low densities to 
primarily collisional at high densities. The low-density limit is ~ 10 cm -3 , while for densities 
higher than 10 4 cm -3 , collisional de-excitation becomes more and more important. Most of 
our sources distribute in a regime where the methods are well applicable. In the low-density 
limit (n e < 10 cm -3 ), the error in density is larger due to the saturation of the relation 
between the line ratio and electron density, so these derived densities should be treated with 
caution. An [S II] ratio of 1.42 corresponds to a density of ~ 10cm~ 3 . Only two objects, 
SDSS J011448.68-002946.1 and SDSS J161809.38+361957.8, have ratios greater than 1.42, 
which places the sources at the low end in the density distribution. 



Notes on individual objects. Among the whole sample, only objects which are above a 
certain S/N ratio in the [S II] line were kept for further analysis (Sect. 2). Among these, we 
then individually re-inspect all objects below a density threshold of n e = 140 cm _3 (logn e = 
2.15 cm -3 , i.e., those in regime "A" in Fig. 5) to check for peculiarities in spectral features, 
and to check the robustness of spectral fitting and thus reliability of the density estimate in 
the low- density regime. 

We comment here on those objects which are extreme in the [S II] ratio, FWHM or 
other spectral features: (1) There is one single broad-line object actually located in the 
'zone of avoidance' in density. This is SDSS J011448.68-002946.1. While its [S II] ratio is 
extreme (R([S II] )= 1.47 ± 0.03), and within the errors beyond the 'low-density' regime, its 
other emission lines are not unusual. The profile of [S II] A6731 deviates somewhat from a 
Gaussian, but there is nothing else very peculiar about it. This is the only outlier in Fig. 5. (2) 
The low-density NLS1 SDSS J092247.03+512038.0 is peculiar in that the peak of the [O III] 
A5007 line is blueshifted with respect to the low-ionization forbidden lines (e.g., [O II], [N II] 

The large velocity 



more than 400 


an s 


Zamanov et al. 


2002) 



A5007 is almost as broad as H/3 b (FWHM 1060 km s -1 and 1250 km s" 1 for [O III] A5007 
and H/?b, respectively), while the low-ionization forbidden lines such as [S II] are narrow 
(FWHMQS II]) = 220 km s" 1 ). It is a strong Fe II emitter with R4570 = 1.32. (3) SDSS 
J161809. 38+361957.8 is the object with the most extreme combination of H/5b line width 
(FWHM(H/?b) = 700 km s -1 for a Lorentzian profile fit) and density (n e = 2.0cm -3 ). The 
emission lines have a very wide range of excitation, the highest corresponding to coronal lines 
of [Fe X] A6374 and [Fe XI] A7892. We also detect the broad low-ionization line O IA8446, 



2000), and that NLR density increases toward small radii (e.g., Fraquelli, Storchi-Bergmann & Binette 2000; 
Bennert et al. 2006a, b). 



-14- 



which is generally produced in a region with very high density (e.g. Komossa & Bade 1999). 
It is the only low-density source of our sample which has the O IA8446 line detected. (4) 
The low-density NLS1 SDSS J120226. 76-012915. 3 shows exceptionally strong Fe II emission 
of R4570 = 2.74. It is the most extreme Fe II emitter in our sample. (5) Finally, we 
mention that among the high-density sources, the spectrum of the BLS1 SDSS J013527.85- 
004448.0 is special in that it shows the highest ratio of [O II] /[O III] of the whole sample 
([OlI]/[OHI] obs = 1.6)@. 

5.2. Density versus FWHM(H/3 b ): a zone of avoidance 

One of our main goals is to examine whether or not there is a difference in electron 
density n e between NLS1 galaxies and BLS1 galaxies, in order to test different NLS1 models. 
The [S II] A6716/A6731 line ratio in our sample ranges from 0.87 to 1.47. We plot the ratio 
versus the FWHM of H/3b in Fig. 5. The typical error on the ratio is ~0.06. The histograms 
of the ratios for the NLS1 and BLS1 samples are also shown. NLS1 galaxies show [S II] 
ratios in the range from 0.94 to 1.43. 17 out of 55 NLS1 galaxies have a ratio higher than 
1.28, while only one out of 39 BLS1 galaxies has a high ratio (i.e., 1.47). The other 38 BLS1 
galaxies occupy the range from 0.87 to 1.27. The average ratios and la deviations of the 
two samples are 1.23 ± 0.12 and 1.12 ± 0.10, respectively. 

In Fig. 5 (lower panel), we display the inferred electron density against the FWHM of 
H/3b. We find that the sources do not homogeneously populate the n e -FWHM(H/3b) diagram. 
The key detection is a 'zone of avoidance 1 in the diagram. While the 38 BLS1 galaxies avoid 
low average densities, and all show n e > 140 cm -3 (regime C), NLS1 galaxies show a larger 
scatter in density in the range n e = 2 ~ 770 cm" 3 , including a significant number of objects 
with low densities. 17 out of the 55 NLS1 galaxies under study show n e < 140 cm -3 (regime A) 
and are clearly separated from the range occupied by BLS1 galaxies. The other 38 NLS1 
galaxies overlap well with the range in density for BLS1 galaxies (regime B)@. The average 
electron densities for zone A, B and C are 69, 294 and 380 cm -3 , respectively. We apply the 
K-S statistical test on the distributions of density for NLS1 galaxies vs BLS1 galaxies. The 
resultant K-S test probability is 0.0002, which means that the two density distributions are 



8 More detailed account on multi-wavelength properties and on (unusual) line profiles and line ratios of 
individual sources will be given in a follow-up paper (Xu et al. 2007, in prep.) 

9 Results are robust, no matter whether a Gaussian or Lorentzian profile is used for the broad component 
of H/3. A larger scatter, rather than a strict cut-off in density for BLS1 galaxies might appear if the sample 
size increases significantly, but we expect our findings do still hold on average. 
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significantly different. 

5.3. [O I] emission versus FWHM(H/5 b ) 

In order to probe also the high-density NLR regime, well above the critical densities 
of the two [S II] lines, we concentrate on the emission line [O I] A6300. Even though 
the ratio [O I] /Ha is also influenced by other parameters, [O I] is strongly boosted for 
high densities (e.g., Fig. 4 and 8 of Komossa & Schulz 1997) and we thus checked for any 
trends and correlations between [O I] intensity and FWHM(H/3b) (Fig. 6), which we regard 
as a supplement to the n e -FWHM(H/3b) diagram (Fig. 5). In Fig. 6, we plot the ratio 
[O I] A6300/Ho; n as a function of the FWHM(H/3b). An anti-correlation between the two 
parameters is found with a correlation coefficient r s = —0.34, and a probability of the null 
correlation P nu ii < 10~ 2 , calculated from the Spearman rank-order correlation analysis. The 
trend meets our expectation from Fig. 5, in the sense that BLS1 galaxies show, on average, 
higher [O I] /Ha than NLS1 galaxies. 

6. Reality of the zone of avoidance 

The sources of our sample do not homogeneously populate the n e -FWHM(H/3b) dia- 
gram, but show a 'zone of avoidance', in the sense that BLS1 galaxies lack low average NLR 
densities. Is there any data analysis or selection effect that could cause a spurious zone of 
avoidance in the BLS1 galaxies regime, or could mimic a larger density scatter in the NLS1 
galaxies regime? We discuss and reject several possibilities in turn. 

Magnitude distribution. Firstly, we note that the magnitude distributions of the NLS1 
and the BLS1 sample are similar (also in dependence of redshift; Fig. 3). We note that both 
extremes of the distribution in the n e -FWHM(H/3 b ) diagram include both quasars and 
Seyferts (but more Seyferts, since the total number of Seyferts is higher than the quasars). 

Profile shape. Secondly, our results are robust, independent of the profile the broad 
emission lines are fit with, either a Lorentzian profile or a combination of two Gaussian 
profiles (see Fig. 5). 

Atmospheric absorption effects. We then checked whether the density estimate for the 
small-FWHM regime is more unreliable (thus, a larger scatter) because at a certain redshift 
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range (z =0.130-0.140), the atmospheric O2 absorption line at 7620A overlaps with one 
or the other sulfur line, making line estimates more unreliable (even though atmospheric 
absorption is generally already corrected for the released SDSS spectra) o In that case, 
sources in the regime with the larger range in densities (zones A & B in Fig. 5) should have a 
specific, narrow redshift range. This is not the case, however. The low density objects vary 
in redshift between z =0.034-0.280, while the rest of the sample shows z =0.038-0.289. 

Signal/Noise. In order to get an accurate [S II] line ratio and thus density measurement, 
we required [S II] AA6716,6731 to have S/N greater than 5 which is quite a common cut-off 
imposed on SDSS data. 90% of our spectra are actually of S/N ([S II]) > 10. In order to see 
whether the S/N of our spectra affects in any kind of way our measurement of the [SII] ratio, 
we have checked whether one class of objects (those with low density) or the other class 
(with high density) shows an excess of sources with low S/N. However, this is not the case. 
We have run a Spearman rank-order correlation analysis on the distribution of the [S II] 
ratio in dependence of S/N, and we do not find any correlation (r s = —0.05, P nu u = 0.63). 

Faint broad wings in the Balmer lines. In case of very broad Balmer lines, if the broad 
wings of the line are over- or underestimated, then a corresponding small under- or over- 
estimate in the [S II] lines might result. The [S II] A6716 line might then be slightly more 
strongly affected, and that might marginally change the measured [S II] ratio. However, in 
most of our sources the [SII] lines are well separated from faint broad wings (see Fig. 2) 

High and low density limits. [S II] A6716/A6731 is only a good density diagnostic in a 
certain range of densities. The low-density limit is ~ 10 cm -3 , while for densities higher than 
10 cm" 3 , collisional de-excitation becomes more and more important. Most of our sources 
distribute in a regime where the methods are applicable. Since the conversion of the [S II] 
ratio into density is non-linear, the estimate of density in the low-density regime (Fig. 5.3 of 
Osterbrock 1985) comes with larger errors. These are shown in Fig. 5. Within these errors, 
objects with small FWHM(H/5b) still have, on average, lower densities. 

The whole diagram (i.e., across the whole FWHM(H/?b) range), may, to a minor extent, 
be biased toward excluding very high densities by the original selection effect of using only 
spectra with significant S/N in the [S II] lines. This pre-selection cannot be avoided. How- 
ever, since [S II] is collisionally de-excited for densities well above 10 4 cm~ 3 , the pre-selection 



10 The other atmospheric O2 absorption line at 6870A will not overlap the sulfur lines, given the redshift 
range 0.034-0.289 of the objects in our sample. 



criterion then is prone to excluding objects with extremely high densities, should they exist 
at all. 



Influence of temperature. The estimate of density using the ratio [S II] A6716/A6731 
also depends on temperature. As described in Sect. 5.1, we fix T e = 10000 K. An estimate 
of the temperature is possible with the temperature-sensitive emission-line ratio [O III] 
A4363/[0 III] A5007. However, [O III] A4363 is a faint line which generally comes with large 
measurement errors. We have inspected the highest-quality spectra available in the [O III] 
A4363 wavelength range (high S/N and/or easy deblending of H7 and [O III]); both, among 
the lowest density object, and, for comparison, a similar number among the higher density 
objects. We do not find systematically higher temperatures for the low-density objects. We 
further note that g iven the scaling of density with temperature, n e (T) = n e (obs) x a/T/10 000 
(jOster brockl Il989h . even an increase in temperature up to 40,000K (close to the average of 
four Syl galaxies of Bennert et al. 2006b), would only change density by a factor of 2, much 
less than the scatter in the observed density values. We come back to this point in Sect. 7.2. 



Starburst contribution. If there was a strong systematic starburst/H II contribution to 
the emission lines, then we would expect the density estimates to be biased toward lower 
values, more typical for H II regions. Since our sample shows unambiguous AGN-like line 
ratios, we expect the vast majority of the emission to come from the NLR, rather than 
star-forming regions. We further use the star-forming indicator, [O II] A3727 (e.g, Hippelein 
et al. 2003; Kewley & Geller 2004; Ho 2005), to track the starburst contribution to the 
optical emission lines of the AGN of our sample. [O II] A3727 is prominent in H II regions, 
while the ratio [O II] /[O III] is observed and predicted to be relatively weaker in Seyfert 
galaxies. We compare the line ratios [O II] A3727/[0 III] A5007 of NLS1 galaxies and BLS1 
galaxies, in order to check whether or not the [O II] emission is on average stronger in 
NLS1 galaxies. In Fig. 7 we show the distributions of the ratio prior to reddening correction 
and after reddening correctional]. No difference is present between the [O II] /[O III] ratio 
for NLS1 galaxies and BLS1 galaxies. Thus, we conclude that the starburst contribution 
is not a possible explanation for the n e -FWHM(H/? b ) distribution. The average observed 
[O H]/[0 III] ratios for the NLS1 and BLS1 sample are 0.32 and 0.29, while the average 
reddening corrected ratios are 0.35 and 0.33, respectively. 



11 For dust reddening correction we use the average reddening curve of Osterbrock (1989, his Table 7.2) 
and an intrinsic value of Ha n /H/3n = 3.1. 
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Spatially resolved NLRs and viewing angle effects. Finally, we checked whether 
nearby objects might be spatially resolved, such that the low-density outer part of the NLR 
is missed in the fiber which would then result in an overestimate of the density in these 
objects. We find that the fiber diameter of 3" corresponds to ~ 2kpc at z = 0.034 (the 
lowest redshift in our sample), while the gas in the NLR is typically distributed over a 
distance r ~ 10 — 1000 pc from the nucleus. Therefore, SDSS fibers would always cover the 
entire NLR. 



7. On the origin of the zone of avoidance 

The key result of this study is the detection of a zone of avoidance in the density - 
FWHM(H/? b ) diagram (Fig. 5): BLS1 galaxies (FWHM(H/3 b ) > 2000 km s -1 ) avoid low av- 
erage densities, and all show n c > 140 cm -3 . On the other hand, NLS1 galaxies show a larger 
scatter in densities in the range n e = 2 ~ 770 cm -3 , including a significant number of objects 
with low densities. The results obtained for [S II] are consistent with the [O I] -FWHM(H/3 b ) 
diagram as shown in Fig. 6, which shows lower average [O I] /Ha in NLS1 galaxies. In the 
following, we first confront our results on density with predictions or indications on density 
effects of existing NLS1 models, and then discuss further possibilities. 



7.1. Coupling between NLR and BLR? 

Given occasional reports of signs of lower or higher-than-average BLR density of NLS1 
galaxies (e.g., Kuraszkiewicz et al. 2000; Rodriguez- Ardila et al. 2000a; Marziani et al. 2001; 
Komossa & Mathur 2001; Xu et al. 2003; Bachev et al. 2004) we then also expect the same 
trend to reflect in the NLR densities - if the properties of NLR and BLR are closely linked^. 
Such models can then be tested by measuring the NLR density. We start with a short review 



12 The exact relation between BLR and NLR in AGN is still under examination. While some models do 
predict a close link, others assume or indicate that both cloud components arc of different and unrelated 
origin. Models that argue against a common link include indications that NLRs are just normal interstellar 
medium (ISM) in the host galaxy, while the BLR has separate origin. Also, the fact that the FWHMs of 
the Balmer lines of BLR and NLR do not generally correlate (e.g., BG92; Grupe et al. 1999; Vaughan et 
al. 2001; Xu et al. 2003) suggests different kinematic components. On the other hand, there are models 
and observations that do predict a link, including common wind outflows (e.g., Schiano et al. 1986); the 
suggestion that the apparent gap between BLR and NLR is solely caused by dust effects (Netzer and Laor 
1993); and observational links between BLR and NLR parameters, in particular, between line strength, line 
asymmetry and the shift of the line centroids (e.g., Xu et al. 2003). 
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of previous measurements, then come back to our new results. 



Few NLS1 galaxies have b een studied with respect to the density of their NLR, so far. 



Rodriguez- Ardila et all (l2000al ) reported NLR density measurements of seven NLS1 galaxies, 
based on the [S II] A6716/[S II] A6731 line ratio. They found lower average density in the [S II] 
emitting zone in NLS1 galaxies than in BLS1 galaxies. Veron-Cetty et al. (2004) presented 
a high S/N optical spectrum of the NLS1 galaxy IZw 1 and concluded that the bulk of the 
NLR is unlike that of most Syl galaxies. It is of unusually low excitation and dominated by 
lines of high critical density, while lines like [O III] A5007 and [S II] AA6716,673 1 are weak, 
hey in fer a density of the low-ionization part of the NLR of n e = 10 6_7 cm~ 3 . lLaor et al. 



(Il997bl ) detected on the same object very weak [C III] A 1907 and [Si III] A 1883 emission in 
a high S/N UV spectrum, suggesting a NLR component with n e ~5x 10 5 cm -3 . 

Rodriguez- Pascual et al. (1997) analyzed the UV properties of a sample of NLS1 galax- 
ies. Based on their photoionization m odeling of the emiss i on lin es, they tentatively favored 



a low-density BLR in NLS1 galaxies. iFerland fe Perssonl (119891 ) presented observations of 



Ca II emission lines from AGN with strong Fe II emission. They suggested that the BLR 
density in strong Fe II emitters is not higher than in other sources, based on the ratio of 
the forbidden lines relative to the Ca II triplet. On the other hand, other studies (e.g., 
Kuraszkiewicz et al. 2000; Wills et al. 2000; Marziani et al. 2001; Bachev et al. 2004) 
favored high-density BLRs, based on the Si III] A1892/C III] A1909 ratio which is sensitive 
to density; such results are consistent with predictions by Gaskell (1985) that UV spectra of 
NLS1 galaxies would show a larger Si III] A1892/C III] A1909 ratio. More recently, in studies 
of the optical-UV emission-line spectra of AGN, Marziani et al. (2001) and Bachev et al. 
(2004) reported indications of a systematic increase in density toward AGN with smaller 
FWHMs of the broad component of H/5 based on the line ratio Si III] A1892/C III] A1909 
(they could not exclude the alternative interpretation of varying metal abundances, though). 
Comastri et al. (1998) reported detection of a deviation of the BLR Balmer line ratio Ha/H/? 
from the recombination value, such that the ratio was below the recombination value. They 
interpreted this as an indication of high density of the BLR of the NLS1 galaxies. 

Our finding is that NLS1 galaxies have lower average NLR density. Therefore, if NLR 
and BLR are closely coupled, then our results favor models which also predict lower average 
BLR density in NLS1 galaxies (i.e., a stronger scatter including objects with lower density). 
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7.2. Supersolar metallicities and temperature effects 

There are several indications that NLS1 galaxies have supersolar metallicities (e.g., 
Mathur 2000; Komossa & Mathur 2001; Nagao et al. 2002; Shemmer and Netzer 2002; 
Romano et al. 2004; Fields et al. 2005). We already noted in Sect. 6 that those objects with 
reliable [O III] A4363 measurements do not indicate enhanced temperatures in the low-density 
objects. However, the number of good spectra is still relatively small. Nagao et al. (2001) 
also examined the temperature-sensitive [O III] ratio for a small sample of NLS1 galaxies 
and concluded that they do not significantly deviate from BLS1 galaxies. On the theoretical 
side, whether an increase in metals, and specific elements in particular, leads to increased 
heating or cooling of the gas, and thus an increase or decrease of its temperature, needs 
to be assessed by detailed photoionization calculations. In general, increasing the oxygen 
abundance first leads to a decrease in temperature since oxygen is an important coolant. 
Such an effect would shift the data points in the n e -FWHM(H/3b) diagram to even lower 
densities (because of the dependence of the [S II] ratio on temperature). In summary, we 
do not expect metallicity effects to play a dominant role in explaining the n e -FWHM(H/3b) 
diagram. 



7.3. Starburst contribution 



If we had a strong starburst contribution in a fraction of our sources, then this would lead 
to lower measured density because HII regions have lower average density (e.g. Osterbrock 
1989; Ho et al. 1997). If NLS1 galaxies are "young objects" (e.g., Mathur 2000; Mathur, 
Kuraszkiewicz & Czerny 2001; Grupe & Mathur 2004; Mathur & Grupe 2005), still in the 
process of growing their black holes, then they may possibly also show enhanced starburst 
activity. However, as already being checked in Sect. 6, NLS1 galaxies, including low-density 
objects, do not show signs of stronger starburst contribution to their optical emission lines 
than BLS1 galaxies. We thus conclude that starburst effects are neglig iblfl 



13 In a study of a sample of 74 post-star bur s t type I AGN, which underwent a strong recent star formation 
epoch but stopped forming stars, Zhou et al. ( 2005^ found that more than half of them are NLS1 galaxies. 
We checked whether their sample, which also includes BLS1 galaxies, shows any special preference for 
certain density values, particularly whether the post-starburst NLS1 galaxies would preferentially have low 
density. We find that their sources populate similar areas in the n -FWHM(H/3b) diagram. We have 
few post-starbursts among the low-density NLS1 galaxies in our sample, based on the lack of evidence of 
strong post-burst features in their host galaxies (equivalent width of H<$ absorption line EW(H<5) > 5 A for 
post-starbursts. 
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7.4. NLR extent 

BLS1 galaxies would lack an observable low- density NLR component, if the low-density 
part of the NLR was selectively obscured. However, spatially-resolved spectroscopy of the 
NLRs of nearby Seyfert galaxies shows that density declines as a function of cloud distance 
from the center. This would imply that it is the outer parts being obscured, leading to a 
peculiar and unlikely geometry. 

Alternatively, lower average NLR density of the small-FWHM(H/3b) objects may imply 
that their NLRs are, on average, more extended, since density declines outward. Indeed, 
there were early suggestions that the emission-line regions of NLS1 galaxies are at larger 
nuclear separations than those of BLS1 galaxies (e.g., Giannuzzo et al. 1999). Such a NLR 
model is consistent with a larger BLR distance in NLS1 galaxies (e.g., Wandel & Boiler 1998, 
Puchnarewicz et al. 2001). 

However, we note, that other parameters also play a role in determining NLR extent. 
In particular, there appears to be a scaling between luminosity of [O III] and NLR extent 
(Bennert et al. 2002, Schmitt et al. 2003), such that more luminous objects have more 
extended emission-line regions. We checked whether our sample shows a correlation between 
density and [O III] luminosity. No obvious trend is found. 

7.5. Fraction of matter-bounded clouds 

A substantial part of the [S II] emission is produced in the partially ionized zone of the 
NLR clouds. If the fraction of matter-bounded clouds varies as a function of the distance 
of the NLR clouds from the nucleus, and if the fraction of matter-bounded clouds in BLS1 
galaxies is higher at larger nuclear distances than it is in NLS1 galaxies, then (given ob- 
servations of density stratification within the NLR from high to low density) some of the 
low-density clouds in BLS1 galaxies partly escape our measurement, with the consequence 
that the average measured density in BLS1 galaxies is higher. Or, to put it the other way, a 
larger scatter in the number of matter-bounded low-density clouds in NLS1 galaxies would 
lead to a larger scatter in their average measured densities. 

There are no direct measurements of the fraction of matter-bounded clouds in NLS1 
galaxies and BLS1 galaxies. However, several observations do indicate a wide of range of 
column densities of the emission line clouds in NLS1 galaxies: Ferland & Persson (1989) need 
very high column densities to reproduce the strength of the Calcium emission in objects with 
strong Fe II emission, while Rodriguez-Pascual et al. (1997) infer a higher fraction of matter- 
bounded BLR clouds in NLS1 galaxies, based on UV observations. Regarding the NLR of 
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NLS1 galaxies, Contini et al. (2003) require matter-bounded clouds in order to reproduce 
the strength of high-ionization iron lines. Photoionization models of the NLR of NLS1 
galaxies presented by Rodriguez-Ardila et al. (2000b) invoke a mixture of matter-bounded 
and ionization-bounded clouds, the inner, high-density clouds mostly matter-bounded and 
producing high-ionization lines; the outer, low-density clouds mostly ionization-bounded and 
producing mostly low-ionization lines. However, Rodriguez-Ardila et al. (2005) discuss a 
different type of NLR models, involving shocks and photoionization, and in that model it is 
shocks which produce the low-ionization lines. 

Future spatially resolved long-slit spectroscopy of the NLRs of nearby NLS1 galaxies 
will allow a direct comparison of the density profiles of BLS1 galaxies (e.g., Bennert et al. 
2006b and references therein) with those of NLS1 galaxies. Making use of other line ratios 
in addition to [S II] may also allow us to determine the fraction of matter-bounded clouds 
as a function of radius, and enable us to measure the NLR extent of NLS1 galaxies. 



7.6. ISM of the host galaxy 



The NLR clouds are most likely directly related to the ISM of the host galaxy. The 
NLR properties then might also reflect different gas enrichment mechanisms of the ISM. For 
instance, the gas could be due to local stellar processes or transported inward from much 
larger scales (or else could arise from the inner parts in case of outflowing winds, driven by 
the central engine; see Sect. 7.7). The absence of low-density clouds in BLS1 galaxies may 
thus reflect the properties of the ISM in the host galaxy. 

The host galaxies of NLS1 galaxies are often spirals, but not much is known about 
their systematic properties. There are indications that the host galaxies of NLS1 galax- 
ies and BLS1 galaxies differ, particularly in the sense that large-scale bar s are more com- 
mon in NLS1 galaxies (jCrenshaw. Kraemer fc Gabelll2003l ; lOhta et al.ll2007l ). and that NLS1 
galaxi es have more grand- design dust s pirals and a higher fraction o f nucle ar star-forming 
rings (jDeo. Crenshaw fc Kraemer! 120061 ) . Moreover, iKrongold et ajj (120011 ) examined the 
host galaxies and found NLS1 galaxies reside in galaxies with smaller diameters than BLS1 
galaxies. Previous studies on mass-luminosity/sigma relation led to conflicting possibilities, 
i.e., NLS1 galaxies are hosted by less luminous galaxies (e.g., Wang & Lu 2001; Botte et 
al. 2004, 2005) or on the contrary, more luminous galaxies for a given black hole mass (e.g., 
Mathur, Kuraszkiewicz & Czerny 2001; Grupe & Mathur 2004; Bian & Zhao 2004; Ryan et 
al. 2007). 



Many of the morphological differences of the host galaxies of NLS1 galaxies and BLS1 
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galax ies are due to the presence or absence of a large-scale stellar bar (IDeo. Crenshaw fe Kraemer 



20061 ). Dynamical models show that a bar potential can efficiently drive gas from the outer 
regions (several kpc) to within ~1 kpc from the nucleus, at which point the bar-driven 
gas flow slo ws or even stalls at the inner Lindblad resonance and the infalling gas will 



form a disk (IShlosman, Bege 



( jSimkin. Su fe Schwarz 



1980; 



man fe Frank! Il990f). Gas inflow a long a galactic stellar bar 



Shlosman. Peletier fe Knapenll2000l ). has been proposed to be 



one possible trigger of AGN activity (e.g. Shlosman, Begeman & Frank 1990). However, 
how these processes affect the ISM density distribution on the scales of the NLR is presently 
unclear. Further studies of ISM properties of NLS1 and BLS1 galaxies on sub-kpc scales are 
needed to further address this issue. 



7.7. Outflows 

Lawrence et al. (1997) suggested that the density of an outflowing wind might be 
an important ingredient in understanding emission-line parameters of AGN and correla- 
tions among them. They loosely suggested strong Fe II emitters have the denser winds, 
in the context of a model where parts of the BLR are mechanically heated and produce 
the Fe II emission. We do find that BLS1 galaxies and NLS1 galaxies differ indeed in the 
density of their NLRs, but that it is actually the BLS1 galaxies which do harbor the higher 
density clouds. A correlation analysis further shows that n e is anti-correlated with R4570 
(r s = —0.47, P nu n < 10~ 4 ; Fig. 8). The low-density objects show larger-than-average R4570 
compared to the high-density NLS1 galaxies. The average R4570 for the low-density objects 
and the high- density NLS1 galaxies is 0.95 and 0.61, respectively. 



Schianol (119861 ) predicted that the average NLR density should be higher in very lumi- 
nous objects than in lower luminosity objects. The prediction was based on their ' quasar 
wind ' model, i.e., the NLR is the result of the interaction of AGN ionizing photons and a 
thermal wind on dense, massive interstellar clouds. However, our result conflicts with this 
prediction. The low-density regime in the n e -FWHM H/5b diagram (Fig. 5) is not dominated 
by low-luminosity systems, but a mixture of high- and low-luminosity NLS1 galaxies. The 
quasars and low-luminosity Syl galaxies have almost the same average NLR density of ~ 
290 cm 



-3 



There are indications that many (but not all; e.g., Xu et al. 2003; Williams, Mathur 
& Pogge 2004) NLS1 galaxies accrete close to or even above the Eddington limit (e.g., 
Boroson & Green 1992; Wang et al. 1996; Boiler et al. 1996; Laor et al. 1997; Sulentic, 
Marziani & Dultzin-Hacyan 2000; Boroson 2003; Grupe 2004; Grupe & Mathur 2004; Collin 
& Kawaguchi 2004). This likely comes with the presence of strong outflows. If these still 
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propagate up into the NLR, then we may expect that the NLR gas in such objects is actually 
more tenuous. Radiation-pressure driven wind models predict a decrease of accretion rate 
with increasing width of the broad component of H/5 (e.g., Nicastro 2000; Witt et al. 1997). 
Among the NLS1 population itself, it should then be the objects with accretion rates closer 
to Eddington that drive the stronger winds and thus have more tenuous, low-density NLR 
components. In order to test this, in Fig. 9, we plot Eddington ratio Lboi/LEdd as a function of 



densi ty n e . We estimate the bolometric luminosities using L bol w 9AL A (5100A) (IKaspi et al. 



20001 ) . while the Eddington luminosities are calculated using the black hole masses deter- 
mined using the BLR radius and the velocity of the BLR gas (e.g. Peterson 1997). We 
find that Lboi/Lsdd is higher in NLS1 galaxies than in BLS1 galaxies. An anti-correlation 
of decreasing electron density with increasing Eddington ratio can be seen across our en- 
tire sample of NLS1 and BLS1 galaxies (r s = —0.42, P nun = 10 -4 ). However, among the 
NLS1 population itself, the low-density objects do not show higher-than-average L bo i/L Edd 
compared to the high-density NLS1 galaxies, which may suggest that higher Lboi/LEdd is a 
necessary but not a sufficient condition to lower density. Independent calculation of L/Lsdd, 
e.g., estimation from X-ray observations, is of great importance to check the trend. 49 out of 
93 objects have X-ray counterparts in the RASS catalogs (e.g. Voges et al. 1999). We had a 
first look at the X-ray data. 40 of them have enough photon counts for an estimation of the 
X-ray slope and thus X-ray luminosity in the ROSAT band. The anti-correlation between 
n e and Lx/Lsdd is even stronger (r s = —0.47, P nu n = 3 x 10~ 3 ), particularly in the sense 
that low-density objects do have higher-than-average L x /L E dd compared to the high-density 
NLS1 galaxies, with average L x /L E dd = 3.2 and 1.2, respectively. Further careful study of 
the X-ray spectra and follow-up X-ray observations of the low-density objects will be crucial 
to understand the cause of the correlation. 

If outflow was a key mechanism to explain the lower average NLR density in NLS1 
galaxies, then we would expect that the density ([S II] ratio) scales with the [O III] outflow 
velocity (blueshift). We checked for both, NLS1 galaxies vs. BLS1 galaxies, and within the 
NLS1 sample (high-density objects vs. low-density objects) whether the density correlates 
with theJO III] (peak) blueshift, and only found a weak correlation (r s = —0.29, P nu u = 
7 x 10~ 3 )o However, it is interesting to note that the peak blueshift of [O III] does strongly 
correlate with L bo i/LEdd ( r s = 0.51, P nu n < 10~ 4 ). This correlation then indicates that 
outflows are more common in objects with high accretion rates. We also checked whether 



14 This does not yet exclude that the whole NLR is in outflow. The [O III] peak blueshift was calculated 
relative to the the low-ionization lines; i.e., [S II] and [N II]. Ideally, one should measure the shift between 
host galaxy absorption lines and NLR emission lines, but most of our spectra are AGN dominated with few 
absorption lines detected. 
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there is a correlation between density and the blueshift of the blue wing of [o since 
a preferred interpretation of blue wings is the existence of outflows (or inflows) combined 
with viewing angle effects (e.g., Boroson 2005). A correlation is seen with r s = —0.43 
(-Pnuii — 1 x 10 -3 ). This correlation shows that outflows are stronger in the low-density 
objects (Fig. 10), even if the bulk of the NLR does not participate in the outflow. 

In summary, we find several indications which point toward a link between NLR density 
and outflows, and we tentatively favor the role of outflows in explaining the difference in the 
NLR density between NLS1 galaxies and BLS1 galaxies. 



8. Summary and conclusions 

We have studied one of the largest homogeneously analyzed sample of NLS1 galaxies in 
order to examine whether or not there is a difference in electron density n e between NLS1 
galaxies and BLS1 galaxies. We employ a powerful diagnostic, the density-sensitive line 
ratio [S II] A6716/A6731, to measure the NLR density. We show that the galaxies do not 
homogeneously populate the n e -FWHM(H/3b) diagram. Our key finding is the detection 
of a 'zone of avoidance' in the n e -FWHM(H/? b ) plane: BLS1 galaxies (FWHM(H/? b ) > 
2000 km s -1 ) avoid low average densities, and all show n c > 140 cm -3 . On the other hand, 
NLS1 galaxies show a larger scatter in densities in the range n e = 2 ~ 770 cm -3 , including a 
significant number of objects with low densities. The results obtained for [S II] are consistent 
with the [O I] -FWHM(H/3b) diagram, which shows higher average [O I] /Ha intensity for 
BLS1 galaxies. 

We investigated a number of different explanations for the 'zone of avoidance' in density. 
We find that supersolar metallicities and temperature effects, a strong starburst contribution 
in NLS1 galaxies, and the effect of NLR extent are unlikely explanations. Consequences of the 
fraction of matter bounded clouds, and different properties of the ISM in the host galaxies, 
can only be further judged with future observations. We find several lines of evidence that 
outflows play a significant role in driving the difference in the NLR between NLS1 galaxies 
and BLS1 galaxies, and favor these as explanation for the zone of avoidance in the density- 
FWHM(Hft) diagram. 
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15 57 out of 93 objects clearly show blue wings. The blueshift of the peak of the blue wing was measured 
against the peak position of the core of the [O III] line. 
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Fig. 1. — Examples of the decomposition of the spectrum into starlight and nuclear contin- 
uum. F\ in units of 10~ 17 erg s _1 cm -2 A -1 is plotted against wavelength in A. In each panel, 
the original spectrum, the power-law continuum of the nucleus, the host galaxy spectrum, 
the Fe II template and the residual spectrum are shown from top to bottom. For clarity, 
the residual spectrum is offset by an additive constant. The three examples are drawn from 
the regimes A, B and C of Fig. 5, defined by the width of the broad component of H/5 and 
electron density. 
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Fig. 2. — Examples of the decompsoition of the H/3+[0 III] and Ha+[N II] emission- line 
profiles. Fa in units of 10~ 17 erg s _1 cm" 2 A -1 is plotted against wavelength in A. The 
narrow emission-lines are fit by Gaussian profiles, while the broad Balmer components are 
fit by either Gaussian (G) or Lorentzian (L) profiles. For Balmer lines the narrow and broad 
components are shown, while for the forbidden lines the resulting narrow line profiles are 
plotted. To illustrate the difference between the data and the fit, the resulting residuals are 
shown at the bottom of each panel. For clarity, the residual spectrum is offset by an additive 
constant. The three examples are drawn from regimes A, B and C of Fig. 5, defined by the 
widths of H/3b and electron density. 
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Fig. 3. — Distributions of i-band absolute magnitude (left) and redshift (right). The open 
histograms plot the NLS1 galaxies, and the shaded histograms plot the BLS1 galaxies of our 
samples. 
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Fig. 4.— Ratio of [0 III] A5007 to the total H/3 flux vs. R4570, the ratio of the flux of the 
Fe II complex between A4434 and A4684 to that of H/? for NLS1 galaxies (filled circles) and 
BLS1 galaxies (open circles). The dashed line and dot-dashed line mark the mean R4570 for 
NLS1 galaxies and BLS1 galaxies, respectively. The large filled circles (pink) represent the 
low- density objects from regime A of Fig. 5. 
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Fig. 5. — Top: [S II] A6716/A6731 intensity ratio vs. FWHM of the broad component of 
H/3 for our sample. Filled and open symbols represent the broad H/5 components modeled 
by Gaussian (G) and Lorentzian (L) profiles, respectively. Squares correspond to QSOs; 
triangles to Seyfert Is. The median error bar is given at the upper left corner. The vertical 
dot-dashed line marks the boundary between NLS1 galaxies and BLS1 galaxies in terms 
of FWHM(H/5 b ). Histograms of the [S II] A6716/A6731 ratio of NLS1 galaxies and BLS1 
galaxies are plotted in the left and right panels, respectively. Bottom: Electron density 
obtained from the [S II] A6716/A6731 ratio in dependence of FWHM of the broad component 
of H/5. Symbols as in the top panel. The arrow points to the location of one outlier which is 
off the plot. The dot-dashed lines distinguish areas populated by: (A) NLS1 galaxies with 
small width of H/? b and low density; (B) NLS1 galaxies small width of H/3b and high density; 
and (C) BLS1 galaxies with large width of H/^ and high density. Median error bars of each 
regime are given. Distributions of the electron density of NLS1 galaxies and BLS1 galaxies 
are plotted in the left and right panels, respectively. 
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Fig. 6. — [0 I]A6300/Ha; n ratio plotted against FWHM of the broad component of H/5. 
Filled triangles are data derived by fitting the broad Balmer components with Gaussian 
profiles, while open triangles are data derived by modeling the broad Balmer components 
with Lorentzian profiles. Large symbols represent the low-density objects from regime A of 
Fig. 5. The [O l]A6300/Ha n ratio is correlated with FWHM(H/? b ) (r s = 0.41, P null < 1(T 3 ). 
The solid line shows the ordinary least-square regression fit to filled triangles, the dotted line 
the fit to open triangles. The dot-dashed line marks the boundary between NLS1 galaxies 
and BLS1 galaxies. 
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Fig. 7. — Histograms showing the [O II] A3727/[0 III] A5007 distribution prior to reddening 
correction (left) and after reddening correction (right). The open histograms refer to NLS1 
galaxies, the shaded histograms to BLS1 galaxies. 
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Fig. 8.— Electron density plotted against R4570 for NLS1 galaxies (filled circles) and 
BLS1 galaxies (open circles). The large filled circles represent the low-density objects from 
regime A of Fig. 5. The sources which are off the plot, are indicated by arrows. The density 
is anti-correlated with R4570 (r s = -0.47, P nuU < 10~ 4 ). 
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Fig. 9. — Electron density plotted against the Eddington ratio, Lboi/LEdd- Filled symbols 
correspond to NLS1 galaxies; open symbols to BLS1 galaxies. Symbols are the same as in 
Fig. 8. A trend of decreasing electron density with increasing Eddington ratio can be seen 
(r s = -0.42, P null = lO" 4 ). 
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Fig. 10. — Electron density plotted against the blueshift of the blue wing of [O III]. Symbols 
are the same as in Fig. 8. 



